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The Harris-Aharony criterion for a statistical model predicts, that if a specific 
heat exponent a > 0, then this model does not exhibit self-averaging. In two- 
dimensional percolation model the index a = — \. It means that, in accordance 
with the Harris-Aharony criterion, the model can exhibit self-averaging properties. 
We study numerically the relative variances Rm and R x for the probability M of a 
site belonging to the "infinite" (maximum) cluster and the mean finite cluster size 
X- It was shown, that two-dimensional site-bound percolation on the square lattice, 
where the bonds play the role of impurity and the sites play the role of the statistical 
ensemble, over which the averaging is performed, exhibits self-averaging properties. 
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I. INTRODUCTION 

The influence of disorder on a phase transition is one of the important problems in the 
theory of phase transitions. In experimental measurements, the thermodynamic proper- 
ties of one (or several) large-sized sample (with respect to the number of molecules) are 
usually studied. Therefore, it is important to know, whether a single large-sized sample 
with quenched realization of impurities can represent the properties of a model? The self- 
averaging properties of the system provide an answer to this question. If some quantity is 
self- averaging, the measurement for a single large sample gives a reasonable value for all 
samples of such size. If the quantity is not self-averaging, an increase in the system size does 
not make the measurement for a single sample representative. 

Let us consider some statistical model with impurities on the ci-dimensional lattice with 
a linear size L, in which we average some quantity X(u) for a certain sample, where u is 
the impurity realization from some impurity ensemble Q. Here we assume that X(u) is the 
exact thermal average for impurity realization u. Let us denote by [X] the average over all 
impurity realizations of ensemble Q: [X] = -P( w )A(u;), where P(u) is the probability 

of impurity realization u>. Then, variance Vx = P{oj)X 2 (uj) — [X] 2 . Let us define 

the relative variance of X: Rx = Txfr- This quantity characterizes the size dependent 
properties of the model. If Rx(L) ~ L x , x = —d the quantity X is said to be strongly 
self- averaging. If x < 0, X said to be weak self-averaging, and if lim Rx(L) — > const, than 

X is non-selfaveraging. We will now discuss self-averaging of a disordered system; therefore 
we examine the averaging over all impurity realizations, instead of self-averaging of a pure 
system, when we study the averaging over the whole statistical ensemble (for example, over 
all possible spin configuration for a spin models). 

The Harris criterion states [Q] that the weak randomness does not change the critical 
behavior of the <i-dimensional second-order phase-transition model if the specific heat index 
a < 0, which corresponds to the pure system correlation length index v > 2/d. As was 
first mentioned by Brout ||, far for critiality where the system size is much greater than 
the correlation length L ^> £, all additive thermodynamical quantities are strongly self- 
averaging. The self-averaging properties of a disordered system near the critiality obeys 
Harris- Aharony (HA) criterion see || and [|J. This criterion states that if randomness is 

irrelevant the system is governed by a pure fixed point and the relative variance Rx ~ L"p , 
where a p and u p are the critical exponents of a pure system. If however, the system is 
governed by a random fixed point, then lim Rx = const. This criterion explains the 

numerical results of many works |5], || |7J. Wiseman and Domany in M show that the 
Ashkin- Teller model with a < is weakly self-averaging. 

Two-dimensional percolation, which can be treated as a g-state Potts model with q — > 1, 
is one of the most popular two-dimensional statistical models with the second-order phase 
transition. The two-dimensional percolation critical index a = — ~ < 0, and, in accordance 
with the HA criterion we can expect that this model exhibits weak self-averaging properties. 
We study numerically the relative variances Rm and R x of the probability M of a site to 
belong the maximum cluster (the analog of magnetization) and the mean finite cluster size x 
(the analog of the magnetic susceptibility). It was shown, that two-dimensional site-bound 
percolation, where bonds play the role of impurities and sites play the role of statistical 
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ensemble over which the averaging is performed, exhibits self-averaging properties. The 
article is arranged as follows. In Section 2 we describe the site-bond percolation model and 
two types of impurities distribution; in Section 3 we discuss the self-averaging criterion and 
its phenomenological derivation; in Section 4 we present the numerical results and Section 5 
is the conclusion. 



II. TWO-DIMENSIONAL SITE-BOND PERCOLATION MODEL AND 
CANONICAL AND GRAND-CANONICAL IMPURITY DISTRIBUTIONS 

Site-bond percolation is a combination of purely site percolation and purely bond perco- 
lation. In purely site percolation, one investigates clusters of occupied sites. In purely bond 
percolation, one investigates clusters of sites, connected by occupied bonds. In site-bond 
percolation one, investigates clusters of occupied sites connected by occupied bonds. Each 
site is occupied with probability p s , and each bond is occupied with probability pb- In critical 
point (p* s ,pl) the correlation length becomes infinite. Yanuka and Engelman || proposed 
approximate formula for a critical curve in (p s , pb) plane 

log(Ps) log(p&) = l 
log(p*) log(pJ) 

The two-dimensional site-bond percolation belongs the same universality class as the two- 
dimensional purely site percolation. Let us define by ric the mean number of clusters of size 
C per lattice site. Then, the probability of a site to belong to the maximum cluster is 

M = n Cmax C max (1) 

and the mean finite cluster size 

E n c C* 

X = E n c C (2) 
Near the critical point (pi, pi) this quantities scales as follows 

M( P s,p b = pt)~(p a -p* s )f ) , p s >p* s} /3 = ± (3) 



and 
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X(Ps,Pb =P* b ) ~ (p s -p* a ) 7 , 7 = 7^ (4) 
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The correlation length £,(p s ,Pb = Pt) ~ (Ps — P* s ) v scaling index is u - .,. 

As mentioned above, each bond is occupied with probability pb and is empty with prob- 
ability 1 — p h ] therefore the total number of occupied bonds on the lattice fluctuates. This 
method of generating impurity configurations is known as Grand Canonical (GC) by analogy 
with Grand Canonical statistical ensemble with a fluctuating number of particles. However 
we can fix the number of occupied bonds Nb and then distribute them randomly on the 
lattice. This method of generating impurity configurations is known as Canonical (C) by 
analogy with the Canonical statistical ensemble with fixing number of particles. 
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III. SELF-AVERAGING CRITERION 

Let us explain the HA criterion via simple phenomenological considerations |J. In this 
section, we use the T temperature as a parameter of the model but keeping in mind, that 
in the case of percolation, the role of this parameter is played by the site concentration p s . 

We characterize every sample (impurity realization) u with size L by a pseudo-critical 
temperature T*(u,L). T*(u,L) fluctuates about its mean value and is averaged over all 
impurity configurations T* = [T*(u,L)]. We introduce the reduced temperature for each 
sample 

L = T -^^ (5) 

In the vicinity of the critical point, the quantity X scales as 

X u (T,L) = ISQ u (i u V*) (6) 

Here, p is the exponent characterizing the behavior of [X] at T* and the thermal scaling 
index y t = - is assumed to be universal for all samples. The form of function is assumed 
to be sample dependent. We assume that (this relation is numerically checked in the next 
section) 

(ST*f ~ (SLY ~ L- d (7) 



So the variance of the argument of the function Q u scales as 



(siL^) 2 ~L- d+ i (8) 



We keep in mind that the scaling relation 0a = 2- ud. So, the variance V x of quantity X 
at the critical point scales as 

V x ~ L 2p (sLL yt ) 2 ~ L 2 "~ d+2 ~ = L 2 " + ^ = L?p + t (9) 



and the relative variance ([X] ~ L p ) is 

Rx = ^ = = L» (10) 

We thus obtain the Harris- Aharony criterion: if the scaling index a = const < 0, then the 
model is self- averaging and the relative variance for quantity X is Rx ~ L~ from simple 
scaling relations. 



IV. NUMERICAL RESULTS 

In our computation, the bond play the role of impurities; therefore, we fix the bond 
concentration pi = 0.875 and consider the site concentration p s as parameter of our model. 



1. we generate (by C and GC methods) the bond configuration u. 
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2. for this bond configuration we generate a set of site configurations (for site we use 
only GC method); For each site configuration we, calculate the probability of a site 
belonging to the maximum cluster and the mean finite cluster size; 

3. We average these quantities over site realizations and obtain mean values M{p s ^p b ) 
and x(Ps,Pb)- We split the set of site configurations into ten series to evaluate the 
numerical inaccuracy AM and A% 

4. we perform steps 1.-3. for another bond configuration. 

The Xu){Ps) dependence for three different bond realizations (u = 1,2,3) is shown in 
Fig.|l]. Let us calculate the pseudocritical site concentration p* s (oo) for each bond realization 
io. We assume, that the mean finite cluster size Xu(Ps) has maximum at the pseudocritical 
point pl{u), and therefore, we approximate the data for Xu near the maximum of Xu> by the 
parabola x^iPs) — a — b(p s — p*(u)) 2 and treat p*{uj) as the pseudocritical concentration for 
bond realization to. The locations of pseudocritical points for bond realizations to — 1,2,3 
is shown in Fig.|l] by vertical lines. 

Let us find the critical point p*. We calculate the mean critical site concentration [p*(£)] 
averaged over 100 bond realizations created by GC and C methods as a function of lattice 
size L. These numerical data and results of approximation 

C- method: [p*(L)} c ~ 0.6519(5) - 0.34(10)L-°- 94 ( 9 ) 
GC - method : \p*{L)) GC ~ 0.6511(4) - 0.39(15)L- L02 ( n ) ^ } 

is shown in Fig|| We take the concentration p* = 0.6515 (averaged by C and GC methods) 
as critical. The variance of pseudocritical concentration behaves as (Sp*) 2 ~ L~ 2 (Fig.||), as 
we assumed in the previous section. 

Let us investigate the self-averaging properties of M and x & t t ne critical point (pi = 
0.875, = 0.6515). 

Now we describe computational procedure that we use to calculate the relative variance. 
Here, we follow the method described |5], ||. First, we note that, for each impurity realization 
uj, instead of the exact value of quantity X w we get some value averaged over site 
configurations with a numerical error 

Here, N s is the number of site configurations (the length of Monte-Carlo run) and t u is 
autocorrelation time. To calculate (5XS) we split the MC sequence of site configurations 
of lenth N s into 10 subsequences and treat each subsequence as independent. We define by 
[. . .] the averaging over impurity - bond configurations. Thus, the error of [X w ] averaged 
over N b bond configurations is 

N) 2 ~E(K1-W 2 ) (13) 
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The total total error ^[Xjj has two contributing terms: from the sample to sample 

fluctuation of exact X u about [X w ] and from the fluctuations of X u , averaged over the finite 
number of spin configurations, about X u for each bond realization u. 



a 



N b N s /r« 



(14) 



Here, iV& is the number of bond configurations. We can calculate 



N b N 3 /r„ 



1 



(sx u y 

by averaging the error (5Xjj over the bond configurations lu G Q. So, we can express the 
relative variance via (5[X tt ,]') and (SX U ) 



R 



V: 



x 



x 



N b (6[X U 



{SX U 



(15) 



The numerical data for relative variance Rm and R x computed by the C and GC methods 
in accordance fllSD , are shown in Fig.[| and |] respectively. We can see that points lies on the 
straight lines in log-log scale. Thus, we approximate the data by power function aL b on the 
interval L e [48, 128]. The results of approximation are plotted in Fig.|5] and |] and placed 
below. As we might expect, for the GC method the relative variance is greater, than for the 
C-method because of the fluctuation of bonds number. 



C : R 



M 



GC : Rm 
C : Ry : 



0.0041(3)L-°- 504 ( 19 ) 

- 0.0147(15)L-°- 50(2) 
i 0.0243(2)L-°' 53(2) 



GC : R x ~ 0.0259(5)L 



-0.51(5) 



We see, that in excellent agreement with HA criterion, the relative variance of measured 
quantities obeys the power low dependence R ~ L~2 = jjz . As a result, we can state that 
this model is self-averaging. 



V. CONCLUSION 



We have shown numerically that the two-dimensional site-bond percolation exhibit self- 
averaging properties at critical point. In our numerical experiments bonds play the role of 
the quenched disorder, and the sites play role of the statistical ensemble. We have found 
find that relative variance scales Rm ~ L~? y R x ~ La, where — | = ^ . We assume that we 
can consider the sites as a quenched disorder and bonds, as a statistical ensemble, and the 
results will be the same. 

We can also expect for the site percolation that, if we "freeze" the state of some sites 
and average over other sites, we get the same weakly self-averaging behavior with respect 
this "frozen" sites, which play the role of quenched disorder. Of cause, the same is valid for 
bond percolation. 
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Figure captions 

Figure 1. The mean finite-cluster size xip) f° r three realization of impurity configurations. 
The critical concentrations p* for each realization cu, are shown by vertical lines. 

Figure 2. The mean pseudocritical site concentration [p*] as a function of lattice size L for 
bond realizations, created by C and GC methods, the results of approximation by a power 
function of L, and the extrapolated critical concentration p* = 0.6515 

Figure 3. The variance {Sp*) 2 of critical concentration as a function of lattice size L for 
C (crosses) and GC (triangles) distributions, and the results of approximation by function 
aL- h . 

Figure 4. The relative variance Km of probability M of a site belonging to the maximum 
cluster as a function of lattice size L for GC (crosses) and C (triangles) distributions and 
the results of approximation by function aL~ b . 

Figure 5. The relative variance R x of the mean finite-cluster size x as a function of lattice 
size L C (crosses) and GC (triangles) distributions and the results of approximation by 
function aL~ b . 
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